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ABSTRACT 

Rest-frame UV spectral lines of star-forming galaxies are systematically offset from the galax- 
ies' systemic redshifts, probably because of large-scale outflows. We calibrate galaxy redshifts 
measured from rest-frame UV lines by utilizing the fact that the mean H I Lya absorption pro- 
files around the galaxies, as seen in spectra of background objects, must be symmetric with 
respect to the true galaxy redshifts if the galaxies are oriented randomly with respect to the 
lines of sight to the background objects. We use 15 bright QSOs at z ss 2.5 — 3 and more than 
600 foreground galaxies with spectroscopic redshifts at z ~ 1.9 — 2.5. All galaxies are within 
2 Mpc proper from the lines of sight to the background QSOs. We find that Lya emission 
and ISM absorption redshifts require systematic shifts of AwL yQ = 



-295I|gkms -1 and 

Auism = 145^35 kms -1 , respectively. Assuming a Gaussian distribution, we put la upper 
limits on possible random redshift offsets of < 220 km s -1 for Lya and < 420 km s -1 for 
ISM redshifts. For the small subset (< 10 percent) of galaxies for which near-IR spectra have 
been obtained, we can compare our results to direct measurements based on rest-frame optical, 
nebular emission lines, which we confirm to mark the systemic redshifts. While our A-uism 
agrees with the direct measurements, our Ai>L ya is significantly smaller. However, when we 
apply our method to the near-IR subsample which is characterized by slightly different se- 
lection effects, the best-fit velocity offset comes into agreement with the direct measurement. 
This confirms the validity of our approach, and implies that no single number appropriately 
describes the whole population of galaxies, in line with the observation that the line offset 
depends on galaxy spectral morphology. This method provides accurate redshift calibrations 
and will enable studies of circumgalactic matter around galaxies for which rest-frame optical 
observations are not available. 
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1 INTRODUCTION 

Two questions that are key to understanding galaxy forma- 
tion and evolution are: How do galaxies get their gas? and 
Where does the gas that is ejected via galactic feedback pro- 
cesses end up? A powerful method to study the interface be- 
tween galaxies and the intergalactic medium (IGM) is to exam- 
ine the spectra of bright background objects for absorption fea- 
tures at the redshifts of foreground galaxies. Different studies have 
looked into the relation between Lya absorbing gas and galax- 
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ies (e.gjLanzetta et al.ll 1995l:IChen et alJll998l:lB~owen et al]2002l: 
I Penton et al J |2002|; I Adelberger e t al.ll2003lT2005bl : ICrighton et al] 
2010; Steid el et al] 120100 . as well as metals and galaxies (e.g., 
Lanzctta & Bowcn 1990; Bergeron & Boisse 1991: Steidel et al. 
1994 1 19971: IChenetalj boOll: lAdelberger et all booj l2005bT : 
Ipieri et al.ll2006l : ISteidel et alfeoioh . A~ crucial requirement for the 
success of such studies is the availability of accurate redshifts for 
the foreground galaxies. To illustrate this point, we note that ob- 
servationally inferred outflow velocities correspond to a change in 
redshift of order Az/(1 + z) = Av/c ~ 10" 3 and that at z = 2 
the difference in Hubble velocity across a proper distance of 1 Mpc, 
which exceeds the virial radius of the typical star-forming galaxy 
at that redshift by more than an order of magnitude, corresponds to 
Az/(1 + z) ~ 7 x 10~ 4 . Clearly, the systematic errors on the red- 
shifts of foreground objects need to be <g 0.1% in order to map the 
physical properties of the gas in and around the haloes of galaxies. 
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The most accurate galaxy redshifts are measured from absorp- 
tion lines arising in stellar atmospheres. Rest-frame optical stellar 
absorption features are routinely detected in spectra of galaxies in 
the local Universe. For high-z galaxies rest-frame UV stellar ab- 
sorption lines are detectable in stacks of ~ 100 high-quality spec- 
tra an d have been used to verify the accura cy of redshift calibra- 
tions ( Shapl ev et al.ll2003l : Isteidel et aljkoioh . One can also mea- 
sure accurate redshifts from nebular emission lines from stellar 
HII regions. Lines such as Ha (A6563), H/3 (A4861), and [OIII] 
(AA4959, 5007) are strong and easily accessible for galaxies in the 
nearby Universe because they are in the rest-frame optical. Measur- 
ing redshifts from these lines is, however, more difficult for high-z 
galaxies. Nebular lines are redshifted into the observed frame near- 
infrared, and at these wavelengths spectroscopy with ground based 
instruments is complicated by the bright night-sky, as well as by 
strong absorption features produced by molecules in the Earth's 
atmosphere. The prospects for near-IR spectroscopy, at least for 
bright objects, will improve when multi-o bject near-IR spectro - 
graphs come online (e.g. Keck I/MOSFIRE; McLean et al.ll2008l) . 
but near-IR spectroscopy will remain challenging for fainter ob- 
jects. 

Alternatively, one may resort to measuring redshifts using 
spectral lines that lie in the rest-frame UV and can thus be observed 
with optical spectrographs. Such lines include the H I Lya emission 
line at 1216A, when present, and UV absorption lines arising in the 
interstellar medium, of which the strongest are C II (A1335), Si IV 
(A1394), and Sill (A1527). However, these lines are not ideal as 
they do not provide accurate redshift measurements. Lya photons 
are resonantly scattered, and the observed line profile depends on 
the kinematics and the amount of gas that is scattering the photons , 
as well as on the dust conten t (e.g.,|V erhamme et al. 2006, 2008; 
Schaerer & Verhammd |200S|: iHa nsen & Oh] 120061 ; iDiikstra et all 



2006al lbl: I Zheng & Miralda-Escudg|2002i) . The observed velocities 



of interstellar absorption lines also depend on the kinematics of 
the gas that produces them. Most high-redshift galaxies observable 
with current facilities have sufficiently high star formation rates to 
drive galactic scale winds. These winds carry some interstellar ma- 
terial out of the galaxies and this causes both the Lya emission line 
and interstellar absorption lines to be systematically offset by a few 
hundred kms -1 from t he galaxy systemic redshifts, as measured 
from nebular lines (e.g..|Pettini et al1l200ll ; [Adelberge r et al.f2003l 
l2005bl ; ISteidel et all2010h~ 

In surveys of high-redshift (z ~ 2 — 3) star- forming galax- 
ies. lAdelberger et aljl2003l l2005bl) and lSteidel etal] b010h dealt 
with these issues by measuring redshifts from rest-frame UV lines, 
obtaining near-IR spectroscopy for a subset (< 10%) of galax- 
ies, and then calibrating the redshifts measured from rest-frame 
UV lines using the rest-frame optical nebular lines. This calibra- 
tion was then applied to the rest of the sample that lacked near- 
IR spectroscopy. However, galaxy populations other than Lyman 
Break Galaxies (LBGs) may require their own calibrations. This is 
a problem as calibration using nebular lines is often not possible. 
For example, most of the galaxies found by future surveys with 
integral field spe ctrographs, such as the Multi Unit Spectroscopic 
Explorer (MUSE. lBacon et al.l2010l) . will only be detected in Lya. 
Also, even for the LBGs at z ~ 2 — 3 it is possible that the subsam- 
ple that is observed in the near-IR is not quite representative of the 
whole sample, for example because galaxies need to be sufficiently 
bright to do near-IR spectroscopy. 

In this paper we present and apply a method for calibrat- 
ing redshifts measured from rest-frame UV lines using the ab- 
sorption features from the surrounding IGM as seen in the spec- 



tra of background QSOs. Several studies have found that the mean 
strength of H I Lya absor ption increases as the distance to the near- 
est gal a xy decreases (e.g.,lchen et alTl 998; Adel berger et alj2003l 
l2005bl. IPieri et alj [20061; iRvan-Webeij l2006t iMorris & lannuzil 
120061 ; IWilman et alj|2007t Isteidel et afl|201Cl) . By using the fact 
that the mean absorption profiles must be symmetric with respect 
to the systemic galaxy redshifts if the galaxies are oriented ran- 
domly with respect to the lines of sight to the background QSOs, 
we determine the optimal average velocity offset for which the ob- 
served mean absorption profile is symmetric around the positions 
of the galaxies along the line of sight (LOS). Although the same 
principle was used to veri f y the accu racy of the calibratio ns pro- 
posed in lAdelberger et al.1 J2003h and Isteidel et al.1 bold) , this is 
the first time that the IGM absorption is used for calibration itself. 
The use of absorption by the IGM makes the method independent 
of the galaxies' luminosities, and does not require rest-frame op- 
tical observations. We find that the required offsets for redshifts 
based on ISM absorption lines and Lya emission lines agree with 
the calibrations inferred by direct comparison with nebular line red- 
shifts. Using this method we also confirm that different galaxy sub- 
samples require different calibrations depending on the galaxies' 
spectral morphology. 

Throughout this work we use f 2 m = 0.258, Q,a = 0.742, 
and h — 0.719 dKomatsu et al.ll2009l) . Unless stated otherwise, all 
distances are proper. 



2 DATA 

The data sample used in this study is p artially (3 out of 15 fields ) 
des cribed in Steidel et a l. (2003; 2004), lAdelberger et alj f2005bl) . 
and Steidel et al. while the full survey will be presented 

elsewhere. Here we give only a brief description. This redshift sur- 
vey was constructed to select galaxies whose U n GR colors are con- 
sistent with star-forming galaxies at redshifts z ~ 1.9 — 2.7. The 
survey was conducted in 15 fields with typical size 5 x 7arcmin, 
centred on bright background QSOs for which there were high- 
resolution spectra suitable for probing Lyman-a absorbing gas 
in the same redshift range. The technique for optically selecting 
galaxies is described in detail in Adelberger et al. (2004) (their 
BX/BM sample). All 659 galaxies used here have been followed 
up spectroscopically with th e optical spectrogra ph Keck I/LRIS- 
B (FWHM « 370 km s" 1 ; Isteidel et alj|2004l) . A subset of 49 
galaxies has been observed with the near-IR spectrograph Keck 
II/NIRSPEC dMcLean et alll 19981) at a higher resolution (FW HM 
« 240 kms -1 ) than achieved by LRIS-B (see lErb et af]|2006l for 
details on this galaxy sample). The sample used here is part of a 
larger sample, with a few thousand galaxies, from which we se- 
lected all galaxies within 2 Mpc from the LOS to the background 
QSOs, and with redshifts that fall in the Lya forest redshift ranges 
of the QSOs in their fields. We define the Lya forest range as the 
part of the spectrum between the quasar's Lya and Ly/3 emission 
lines, excluding the region < 5000 kms" 1 from the quasar's red- 
shift (to avoid the proximity zone). 

The QSO observations were conducted between 1996 and 
2009 with the Keck I/HIRES echelle spectrograph. The spectra 
have a resolution of Av w 7.5 kms - (R ~ 40,000) and were 
rebinned into pixels of 0.04 A. They cover the wavelength range 
~ 3100 — 6000 A. The reduction was done using T. Barlow's MA- 
KEE package, and the spectra were normalized using low order 
spline fits. The typical S/N in the Lya forest region lies between 40 
and 160. 
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Repeated observations of the same galaxies suggest that the 
typical statistical measurement uncertainties are ~ 100 kms -1 
for redshifts measured from ISM absorption lines (zism), ~ 
50 km s - 1 for redshifts measured from Lya emission lines (zLya), 
and ~ 60 kms for redshifts measured from nebular emission 
lines (Zncb)- The systematic offsets between these redshifts and the 
true (i.e. systemic) redshifts, which is the subject of this paper, may 
of course be larger. 



3 METHOD 

Absorption spectra of background objects (QSOs, GRBs, or star- 
forming galaxies) can be used to study the intervening IGM. Ob- 
servations generally show an increase in the mean Lya and metal 
absorption as the distance to the nearest galaxy decreases. One ex- 
pects to find an average absorption profile that is symmetric around 
the galaxies' positions along the LOS. This is true even if the IGM 
is not distributed isotropically around galaxies, which can for ex- 
ample occur if the mean absorption at a given distance from a 
galaxy correlates with the galaxy's inclination. Assuming that the 
galaxies are oriented randomly with respect to the LOS, we still 
expect to see a symmetric absorption profile when averaged over 
an ensemble of galaxiefl If the measured redshifts of galaxies are 
systematically offset from their systemic redshifts, then the mean 
profile will be symmetric around some other velocity point, offset 
from the observed redshift by an amount equal to the systematic 
shift in galaxy redshifts (see Fig.[TJfor an illustration). In this paper 
we use this principle to calibrate the systematic offsets between the 
observed and systemic (i.e. true) redshifts for the galaxies in our 
sample. 

One could imagine that instead of measuring the offset from 
the average absorption profile, we could measure offsets for indi- 
vidual galaxies by observing the IGM absorption features around 
them, and then find the average of such estimated offsets. How- 
ever, near the redshift of a galaxy one often finds several lines at 
different velocities, of which some are associated with the galaxy in 
question, but the others are produced by, for example, a cloud that 
is associated with a neighboring galaxy or an intervening cloud. If 
the absorption line associated with such a cloud were stronger than 
that due to the gas associated with the galaxy, then we would in- 
fer the wrong offset. Also, at small impact parameters even a cloud 
that is associated with the galaxy in question could potentially have 
kinematics that is not dominated by gravity (e.g., it could be a part 
of the galactic outflow); and at any point in the forest there could 
be absorption completely unrelated to the observed redshift. This is 
not a problem when using the average absorption profile, given that 
it is equally likely for such a line to be at a positive as at a negative 
velocity with respect to the galaxy. 

We use Lya rather than metal lines, because the metal-line 
absorption is usually only enhanced within a couple of hundred 
kpc from galaxies, and our sample does not include a sufficiently 
large number of galaxies with such small impact parameters (see 
the discussion below on the required number of galaxies). However, 

1 As the blue side corresponds to a lower Lya redshift than the red side, 
and the mean Lya absorption at a random location is known to increase with 
redshift, we do expect the intergalactic Lya absorption to be slightly weaker 
on the blue side. However, the magnitude of this effect (th e mean trans- 
missi on varies by 0.002 over 10 3 kms -1 at z = 2.3; e.g. Schaye et al. 
120031) is negligible compared to the enhancement in the mean absorption 
near galaxies. 
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Figure 1. The mean flux profile is symmetric around galaxy positions if 
the observed redshifts (z oos , vertical red line) are equal to the true redshifts 
(left panel). If the observed redshifts are systematically offset from the true 
redshifts, the mean profile will be symmetric around some other redshift 
point, offset from the observed redshifts (right panel). 

at higher redshifts, where the Lya forest starts to saturate, it might 
be advantageous to use metal lines, as long as the survey contains 
sufficient galaxies. 

We compare the "blue" and "red" sides of the median Lya 
absorption profile, within ±500 kms -1 of the galaxy redshifts 
of all the foreground galaxies with impact parameters, b, smaller 
than 2 Mpc. We try a large number of velocity offsets, z bs — > 
Zobs ± — (1 + z oos ), where z a bs is the observed redshift, and c 
is the speed of light. Note that a positive Av indicates that the ob- 
served redshifts systematically underestimate the true redshifts.We 
verified that there is no justification for using a redshift dependent 
offset - either because our redshift range is too small for evolution 
to matter (z ~ 1.9 — 2.6), or because there is no dependence on 
redshift in this galaxy population. 

The analysis involves the following steps: i) we assume that 
the observed redshifts are systematically offset from the systemic 
redshifts by An and assign each galaxy a new redshift z' ohs = 
Zobs + — (1 + z bs); ii) for each galaxy we shift the QSO spectrum 
into the galaxy's rest-frame; Hi) for each galaxy we interpolate the 
flux in the QSO spectrum within ±500 km s -1 from z' ohs in ve- 
locity bins of 10 kms -1 to get an absorption profile; iv) we find 
the mediar0 flux profile as a function of velocity separation from 
galaxies; v) we compare the "blue" (—500 — Okms -1 ) and the 
"red" (0 — 500 kms -1 ) side of the median flux profile using \ 2 
statistics: 

2 y (Bi - Ri) 2 
X ti^ 2 (Bi)±^(R0 

where and are the blue and red median flux arrays, a (Bi) and 
er(Ri) are the errors on the median flux, and the sum extends over 
all velocity bins within ±500 kms -1 from the assumed redshifts. 
The errors on the median flux profiles are estimated by bootstrap 
resampling. We recreated the galaxy sample 1000 times by divid- 
ing each QSO field into redshift bins of 1000 kms -1 and then ran- 
domly drawing redshift intervals to recreate each field, allowing 
for individual intervals to be selected multiple times. Each boot- 
strapped galaxy sample consists of galaxies whose redshifts fall in 
the selected redshift intervals. We chose an interval of 1000 km s~ 
because this length is la rger than the correlation length of LBGs 
dAdelberger et alj|2005al) . These five steps are repeated for a large 
number of assumed redshift offsets, yielding a reduced \ 2 curve 
(see Fig. |2j. We expect the \ 2 t° reach a (local) minimum for the 
offset that brings the measured redshifts closest to the true redshifts. 
Note that the x 2 will also have a small value if the applied redshift 
offsets are very large, because in that case the absorption will no 

2 We verified that using the mean rather than the median gives nearly iden- 
tical results for samples with enough galaxies, but that using the median is 
more robust for small samples. 
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longer be correlated with the assumed locations of galaxies, the me- 
dian flux profile will be flat, and the blue and red sides will agree. 
We therefore search for \ 2 minima only between the inferred x 2 
maxima (see Fig. [2] for the characteristic shape of the \ 2 curves). 
In practice we do this as follows: i) we smooth the \ 2 curve with 
a boxcar average with 50 kms -1 width; ii) we find the global \ 2 
maximum and determine the FWHM of the \ 2 peak it belongs to; 
Hi) we search for the second (local) \ 2 maximum that is at least 
one FWHM away from the first one; iv) we find the \ 2 minimum 
between these two peaks, in the un-smoothed curve. 

The value that we chose for the maximum impact parameter, 
b = 2 Mpc, is a compromise between the need for a strong ab- 
sorption signal (the mean absorption is stronger for galaxies with 
smaller impact parameters), and the need for a large number of 
galaxies (there are more galaxies with larger impact parameters). 
We verified that using smaller maximum impact parameters gives 
consistent results. Any maximum impact parameter in the rang^fl 
0.5-2 Mpc yields redshift calibrations that agree with those ob- 
tained for the 6 < 2 Mpc sample within the estimated errors. How- 
ever, the required number of galaxies does depend on the adopted 
maximum impact parameter. For b < 1 and b < 2 Mpc, we obtain 
converged results if we use > 100 and > 200 galaxies, respec- 
tively, although the error bar on the inferred systematic shift con- 
tinues to decrease slowly if more galaxies are included. We also 
varied the velocity range around each galaxy. In this case there 
is a tradeoff between the strength of the signal and the number 
of pixels available for statistics. We found that the velocity range 
of ±500 kms" 1 generally gives the tightest confidence intervals 
around the estimated offsets, but the results are insensitive to the 
precise velocity range chosen. 

We apply the above method to our sample and estimate sys- 
tematic velocity shifts for the subsamples listed in Fig. [3] We make 
a distinction between the whole sample and the NIRSPEC subsam- 
ple due to the poss ibility that this s ubset is not representative of the 
whole sample (see lErb et alj|2006l for details on the sample selec- 
tion and possible biases). 



4 RESULTS 

In Fig.f2]we plot the reduced \ 2 as a function of the assumed veloc- 
ity offset. We find that the best offset for redshifts measured from 
Lya emission is A«Lya = — 295I35 kms -1 (i.e., the Lya emis- 
sion line is systematically redshifted with respect to the systemic 
velocity), while for redshifts measured from ISM absorption lines 
we find Avtsm = I45I35 kms -1 (i.e., the ISM absorption lines 
are systematically blueshifted). 

We estimated the confidence intervals around the best esti- 
mates for A«ism or Ai>L ya using the bootstrap method, described 
in Section [3] We perform all the steps described in Section [3] on 
samples created by bootstrapping the original sample, resulting in a 
best-fit An for each bootstrap realization of the galaxy sample. The 
error on the estimate of Av is the la confidence interval around the 
median Av from the 1000 bootstrap realizations. These la confi- 
dence intervals are shown as the light blue shaded regions in Fig. [2] 
together with the \ 2 curve. 

We prefer to determine confidence intervals by using bootstrap 
resampling rather than A% 2 because the former is more robust. In- 
deed, we found that for small galaxy samples the \ 2 curves become 

3 There are insufficient galaxies with b < 0.5 Mpc to test smaller maxi- 
mum impact parameters. 



too noisy for A% 2 estimates to work, while for large samples they 
yield errors slightly larger than those estimated from bootstrap re- 
sampling. Contrary to A\ 2 , bootstrap resampling does not require 
the errors on the median flux profile to be Gaussian. It also does not 
require the different velocity bins to provide independent measure- 
ments, which they in fact do not because the individual absorption 
lines are broader than the bins in Av. Bootstrap resampling merely 
requires that the redshift regions into which we divide the fields are 
independent, which is true since, as already mentioned, they are 
larger than the galaxies' correlation length. We thus only use the 
X 2 curve to find the Av that minimizes the difference between the 
blue and red sides. 

For the subset of galaxies with near-IR data, the nebular lines 
provide accurate markers of the systemic redshifts of the individ- 
ual galaxies. Comparing these with the rest-frame UV lines, we 
find median offsets of Av-L ya = — 4O6I3Q kms" 1 and A«ism = 
166^35 kms" 1 (based on 42 and 86 galaxies, respectively, includ- 
ing objects that are not within 2 Mpc of the LOS to a QSO). 
The quoted uncertainties correspond to the la confidence interval 
around the median, estimated by bootstrap resampling the galaxy 
sample 1000 times. These results ar e shown as r e d poly gons in 
Fig. [2] For reference, we note that ISteidel et alj bold) quoufl 
mean offsets of Av^ ya = —445 ± 27 kms" 1 and Auism = 
164±16kms -1 . While our result for Avism agrees with the direct 
estimate, the best estimates for AnL yQ differ significantly. 

We therefore employ our IGM calibration method for the 
subsample with nebular redshifts. For Lya the resulting offset is 
Avhya = — 455tg5 kms -1 . This is in good agreement with the 
direct estimate (AvLya = —406^30 kms -1 ) from nebular lines 
for the same galaxies. Both disagree with our measurement for the 
full sample (Afj LyQ = -295±35 kms -1 ). 

We emphasize that our measurement of Avhya for the near-IR 
sample is based on only 26 galaxies, which is generally insufficient 
to obtain a converged result. Indeed, we find that our method fails 
for the subsample of 48 galaxies with both ISM absorption and 
nebular emission lines. The reduced x 2 curve is noisy and does 
not show two clear peaks, which causes the best-fit velocity offset 
to vary strongly with the velocity interval over which we measure 
the median flux. Although the x curve for the subsample of near- 
IR galaxies with Lya emission does look normal and we found 
the result to be insensitive to the precise velocity interval used, the 
number of galaxies may well be too small to obtain a robust mea- 
surement. 

Galaxies whose redshifts have been measured only from their 
Lya emission feature tend to be those with poorer continuum S/N, 
due either to faintness or to poor observing conditions. Such galax- 
ies are under-represented in the near-IR subsample because their 
redshifts are considered less secure and t herefore were s elected 
against for near-IR spectroscopic follow-up terb et alT20 06). It has 
also been observed that the velocity offset o f Lya emission is an ti- 
correlated with the Lya equivalent width ( Shap lev et al.1 120031) - 
stronger Lya lines are expected to have velocities closer to the 
galaxy systemic redshift. These two effects are likely responsible 
for the differences between the full sample and the smaller near-IR 
subsample. 

This difference sends a warning that describing a galaxy pop- 
ulation with a single number is not the best strategy given that the 



4 ISteidel et all <2010l) also quote Ati Lya = -485 ± 175 kms" 1 and 
Ai>lSM = 166 ± 125 kms -1 , but these errors are standard deviations 
around the best offset estimates. 
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Figure 2. The reduced \ 2 obtained by comparing the red and blue sides of the median Lya absorption profile, as a function of assumed velocity offset 
Au = c(z ga i — £ bs) /(I + ^obs). where z ga i and z b s are the true systemic and observed redshifts, respectively. The red and the blue sides are compared 
within ±500 km s~ 1 of all galaxies with impact parameters smaller than 2 Mpc. The left panel shows the results for redshifts measured from Lya emission 
lines (321 galaxies), and the right panel for redshifts measured from ISM absorption lines (590 galaxies). The dashed vertical line shows the (local) minimum of 
the x 2 curve. The vertical, light blue regions show the 1 a confidence intervals, as estimated by bootstrap resampling the galaxy samples. The \ 2 becomes small 
for very large offsets because in that case we are comparing "no signal with no signal" since we completely miss the galaxies. The horizontal red polygons 
show the la confidence interval for the median offset for the near-IR subsample, obtaine d by direct compari son with nebular lines (see al so lSteidel et alj 
l2010l) . The grey and dark blue horizontal lines show the estimates for the mean offset from Steid el et alj feOld) and lAdelberger et alj (2005b) obtained from 
nebular lines, respectively, with the size of the bars indicating the lcr scatter among the galaxies. Our measurements agree well with direct estimates based on 
comparison with nebular lines for ISM absorption lines, but for Lya emission lines they differ significantly. 



line offsets depend on galaxy spectral morphology. Thus, conclu- 
sions drawn from the near-IR subsample may not be generally ap- 
plicable and separate redshift calibrations are needed for each sam- 
ple of interest. 

However, we do confirm the finding from lSteidel et all feoid) 
that absorption line offsets are less sensitive to galaxy spectral mor- 
phology thanLy_Q_£mi^on_Tme offsets. Th is justifies the strate gy 
taken in lAdelberger et alj J2003ll2005bl) and lSteidel etafl J2010h to 
use redshifts measured from ISM absorption lines when they are 
available (~ 90% of the sample). 

As an additional test of our IGM calibration method, we mea- 
sure the systematic velocity offset for redshifts measured from neb- 
ular emission lines, which are known to be close to systemic posi- 
tions. The resulting offset, which is based on a sample of 49 galax- 
ies, is An nc buiar = — 75^30 kms -1 , i.e. consistent with zero, as 
expected (Fig. [3j. 



4.1 Random offsets 

We have demonstrated that one can use the fact that Lya absorption 
by the IGM is correlated with proximity to galaxies to calibrate the 
redshifts inferred from Lya emission and ISM absorption lines. In 
other words, we measured the systematic offsets in the observed 
redshifts. We will now use the same idea to constrain the statis- 
tical (i.e. random) offsets. Statistical scatter in individual galaxy 
redshift offsets may be caused by measurement and instrumental 
errors, but in our case it is probably dominated by the variations 
in the intrinsic galaxy properties. By comparing redshifts inferred 



from rest-frame UV lines with those from nebular lines for the sub- 
set of 89 galaxies that have been observed in the near-IR, one not 
only finds a large systematic offset, but there is also a significant 
scatter of « 180 km s -1 for Lya redshifts, and w 160 km s -1 for 
ISM redshifts. 

We use all galaxies with impact parameters smaller than 2 
Mpc, and measure the median absorption around them. We then 
repeat the procedure after applying a random redshift offset to the 
redshift of each galaxy, drawn from a Gaussian distribution with 
mean zero and standard deviation a. We do this many times for in- 
creasing values of a. We then compare the resulting median flux 
profiles with the original (i.e. a = 0) median absorption profile. 
We expect that the two profiles will be consistent with each other as 
long as the added random offsets are small compared with the true 
random offsets. Using this method we obtain lcr upper limits on 
the random redshift offsets of < 220 km s -1 (< 300 km s -1 ) for 
Lya redshifts, using mean (median) statistics, and < 420 km s -1 
(430 km s -1 ) for ISM redshifts, which is consistent with the in- 
dependent estimates based on the subsample with redshifts from 
nebular lines. The quoted values are upper limits because we mea- 
sure the amount of scatter needed to cause a significant difference 
with respect to the original absorption profile, i.e. the value we find 
exceeds the actual scatter that is present in the sample. 



5 SUMMARY & CONCLUSIONS 

Measuring accurate redshifts for high-jz galaxies is a daunting task. 
Distant objects are faint, making even ground-based rest-frame UV 
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Figure 3. Different subsamples (first column) together with the number of galaxies in the sample (second column) and the best-fit velocity offsets, Av, and 
lex confidence intervals (third column). These best offsets (black diamonds) together with the 1<t confidence intervals (grey horizontal stripes) are shown 
graphically on the right. The results for the full samples are shown as vertical gray lines. 



observations challenging. Added to this is the complexity of the 
gas kinematics and radiation transport, which may cause absorp- 
tion and emission features to be systematically offset from the true 
systemic redshifts. Here we have used a sample of 15 absorption 
spectra of bright background QSOs, with more than 600 foreground 
galaxies within 2 Mpc from the LOS to these QSOs, to probe the 
IGM close to galaxies at z ~ 2 — 2.5. All galaxies have spectro- 
scopic redshifts measured from either the HI Lya emission line 
and/or rest-frame UV ISM absorption lines. In addition, for a sub- 
set of 49 galaxies redshifts have been measured from rest-frame 
optical nebular emission lines, which are thought to be unbiased 
with respect to the systemic redshifts. By utilizing the fact that the 
mean absorption profiles must be symmetric with respect to the true 
galaxy redshifts if the galaxies are oriented randomly with respect 
to the LOS, we calibrated the redshifts measured from rest-frame 
UV spectral lines. 

Our results are summarized in Fig. [3] We found that observed 
Lya emission redshifts require a systematic shift of AuL yQ = 
— 295I35 kms -1 (i.e. the Lya lines are redshifted with respect 
to the systemic velocity) and that ISM absorption redshifts need 
an offset of A«ism = 145l35kms _1 (i.e. the ISM lines 
are blueshifted). For the nebular lines we found Au nc b u iar = 
— 751|q kms -1 , which is consistent with zero, as expected. 

We compared our calibrations to those obtained by direct 
comparison with nebular emission lines, which should be close 
to systemic, for a subsample observed with NIRSPEC/Keck. 
While our measurements for Auism agree, the direct measure- 
ment for the systematic offset of the Lya emission line, A«L ya = 
—406^30 kms -1 , differs significantly. However, when we applied 
our method to the NIRSPEC subsample we found AuLya = 
— 455tg5 kms -1 , in agreement with the direct measurement. The 
discrepancy between the Lya offsets between the near-IR subsam- 
ple and the full sample reflects the fact that the line offset depends 
on galaxy spectral morphology, which is different for the NIRSPEC 
subsample due to the way in which these galaxies were selected. 

After having demonstrated that one can use IGM absorption 
to measure systematic galaxy redshift offsets, we also estimated 



upper limits on the random offsets for individual galaxies. In our 
case these random offsets are likely dominated by intrinsic scat- 
ter rather than by measurement errors. We applied random veloc- 
ity shifts drawn from a Gaussian distribution of varying width to 
the redshift of each galaxy until the median Lya absorption pro- 
file close to galaxies was significantly affected. This procedure 
yielded la upper limits on the random offsets of < 220 kms -1 
(< 300 km s -1 ) for Lya redshifts, using mean (median) statistics, 
and < 420 kms -1 « 430 kms -1 ) for ISM redshifts, which are 
consistent with direct measurements for the subsample with nebu- 
lar redshifts. 

These results will for example be of interest for future Lya 
emitter surveys for which follow-up near-IR spectroscopy will not 
be available. Such surveys would otherwise have to rely on pre- 
vious calibrations for LBGs at z ~ 2 — 3, which might not be 
appropriate for different galaxy populations or redshifts. Without 
applying our self-calibration technique, the redshifts based on Lya 
would likely be systematically wrong by a few hundred kms -1 , 
which would make it very difficult to study the galaxies' environ- 
ments through absorption line spectroscopy of nearby background 
QSOs or GRBs. 

In future papers we will also use the calibrations presented 
here to study the IGM near the galaxies in our sample. 
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